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ABSTRACT 

This work investigates the effects of compressibility on a stationary mode of instability of 
the three-dimensional bound- iy layer due to a rotating disc. The aim is to determine whether 
this mode will be important in the finite amplitude destabilization of the boundary layer. 
This stationary mode is characterized by the effective velocity profile having zero shear stress 
at the wall. Triple-deck solutions are presented for an adiabatic wall and an isothermal wall. 
It is found that this stationary mode is only possible over a finite range of Mach numbers. 
Asymptotic solutions are obtained which describe the structure of the wavenumber and the 
orientation of these modes as functions of the local Mach number. The effects of nonlinearity 
are investigated allowing the finite amplitude growth of a disturbance close to the neutral 
location to be described. The results are compared with the incompressible results of P. Hall 
(Proc. R. Soc. Lond. A406, 93-106 (1986)) and S. 0. MacKerrell (Proc. R. Soc. Lond. 
A413, 497-513 (1987)). 
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1. INTRODUCTION 


The problem considered here is that of the instability of the compressible boundary 
! layer on a rotating disc. The flow due to a rotating disc has a three-dimensional boundary 

layer and exhibits a crossflow instability. This boundary layer is similar to that which 
occurs on swept wings and so can be used to predict the possible modes of instability 
of boundary layers on swept wings. Interest in these types of instabilities has increased 
recently due to the development of laminar flow wings. Our aim here is to investigate 
the effects of compressibility on a mode of instability which will be described below. For 
incompressible three-dimensional boundary layers this particular mode of instability may 
play a crucial role in the finite amplitude destabilization of the boundary layer. 

| One of the first experimental and theoretical investigations of instabilities of three- 

dimensional boundary layers was conducted by Gregory, Stuart and Walker (1955). Using 
the china clay visualization technique they observed a set of stationary vortices in the form 
of spiral streaks rotating with the disc. The number of vortices observed was 28-31 with 
the angle between the tangent at a given point on the vortex and the normal to the radius 
vector being 13°. Linear stability analysis was used by Gregory, Stuart and Walker (1955) 
to explain the nature of the instability. This was described as an inviscid disturbance 
with an effective velocity profile having zero velocity at an inflexion point. The theory 
presented enabled the orientation and the wavenumber of the vortices to be predicted. 
The orientation of the vortices predicted by this theory was found to be consistent with 
the experimental observations but the vortex wavenumber predicted was about four times 
the value observed. 

Later linear stability investigations following on from that of Gregory, Stuart and 
Walker (1955) of the incompressible flow due to a rotating disc were given by Cebeci and 
Stewartson (1980), Malik, Wilkinson and Orszag (1981), Mack (1985) and Malik (1986). 
The results of Malik, Wilkinson and Orszag (1981) show that the effects of streamline 
curvature and Coriolis force may be important. (This result was first found in the context 
of the Ekman boundary layer, by Faller and Kaylor (1966) and Lilly (1966).) Malik (1986) 
obtained the finite Reynolds number continuation of the inviscid mode of Gregory, Stuart 
; and Walker (1955) and showed that at large Reynolds numbers an additional stationary 

' mode can occur. The structure of this mode was subsequently investigated by Hall (1986) 

* who showed that it has a triple-deck structure and that it corresponds to an effective ve- 

locity profile having zero shear stress at the wall. Later MacKerrell (1987) investigated 
the nonlinear development of this stationary mode and obtained the important result that 
it is destabilised by nonlinearity and thus may cause the premature growth of instability 
i waves if it is forced above a certain threshold level. 

i 
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Visible evidence that an additional stationary mode of instability may occur was given 
by the experiments of Fedorov et al. (1976). They observed a stationary vortex structure 
having 14-16 vortices with the angle between the tangent on a point of the spiral and the 
normal to the radius vector being 20°. Thus this structure has a much lower wavenumber 
than that of the Gregory, Stuart and Walker mode and also a different wave angle. 

The main aim of the present paper is to see how the stationary mode of Hall (1986) 
develops into the compressible regime. Most importantly we want to see how far into the 
compressible regime the potentially important subcritical nature of the instability persists. 
We show that the mode of instability is in fact possible only over a finite range of Mach 
numbers. 

In Section 2, we give the full compressible equations governing the flow of interest and 
the boundary conditions to be satisfied. These equations are solved in Section 3 for large 
Reynolds number using a Dorodnitsyn-Howarth transformation to give the basic three- 
dimensional compressible flow due to a rotating disc. We consider two cases of an adiabatic 
disc and an isothermal disc. In Section 4, this basic flow is perturbed by adding a small 
three-dimensional disturbance to it and a nonlinear stability analysis is presented. The 
triple-deck structure of the disturbed flow is described. Although the triple-deck scalings 
here differ from those of the classic triple-deck the same principles apply. A solvability 
condition on the linear lower deck problem produces an eigenrelation, the solutions of 
which are discussed. In order to obtain an equation for the disturbance amplitude we 
perturb the solution a small distance from the neutral location. A solvability condition in 
the lower deck problem produces an equation for the disturbance amplitude, the solution 
of which is discussed. Finally, in Section 5 the effects of compressibility are discussed and 
the results are compared to those obtained for an incompressible flow by Hall (1986) and 
MacKerrell (1987). 

2. THE GOVERNING EQUATIONS 

We will consider the steady flow of a compressible viscous fluid above an infinite flat 
disc rotating with an angular velocity Cl about the 2 -axis. Thus, the fluid motion is induced 
by the rotation of the disc and the effects of streamline curvature and Coriolis force must 
be taken into account. The problem is formulated in terms of cylindrical polar coordinates 
to take advantage of the axisymmetric motion which occurs in the basic flow. Thus, we 
introduce the non-dimensional coordinates and velocities ( r,0,z ) and (v r , vg, v z ), with the 
2 -axis being normal to the surface of the disc. The lengths r and z and the velocities have 
been made dimensionless with respect to some reference length £ and Cl£, respectively. 

Since the fluid is compressible the variation of the density of the fluid must be consid- 
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ered. Hence, the Navier-Stokes equations and the continuity equation are not sufficient to 
solve the problem. Thus, we must also have an equation of state and an energy equation. 
The Reynold number for the flow is defined by 

MX 


R = 


Mo 


( 2 . 1 ) 


where p^ and p^ are the free-stream values of the density and the shear viscosity of the 
fluid, respectively. The Reynolds number is taken to be large in the following analysis. 
The density and temperature of the fluid are non-dimensionalized with respect to their free 
stream values, while the shear viscosity and the bulk viscosity are both non-dimensionalized 
with respect to the free-stream value of the shear viscosity. Finally, the pressure is non- 
dimensionalized with respect to p^ftft 2 . The continuity and Navier-Stokes equation for a 
steady compressible fluid in a reference frame rotating with angular velocity ft are 
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(2.2a, b,c,d) 


Here p is the density, p is the fluid pressure, p and p are the shear and bulk viscosities, 
respectively and V • v is defined by 


Id . ldvg dv z 

v. 2 = 7 -K) + ~ + — . 

We consider the fluid to be a perfect gas then the equation of state is 

1 


P = 


iMl 


pT , 


(2.2e) 


(2.3a) 
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where T denotes the temperature of the fluid. The ratio of specific heats of the fluid is 
defined by 

7 = (2-36) 

c v 

where c p is the specific heat at constant pressure and c v is the specific heat at constant 
volume. The free-stream Mach number is defined by 

M 


M a 




(2.3c) 


where Too is the free-stream temperature and k is the gas constant, which can be expressed 
as k — c p - c„. Using (2.3) the energy equation for a steady flow may be written in terms 
of the temperature and is given by 


(7 - 1 Wi, 


( dT vgdT dT\ dp v 9 dp dp $ 

("'aT + 7 Te + ' v 57j = a? + 7ai + v ‘ai + r 


1 .1 d . , dT. Id., dT. d ., dT., 

' R{i - l)Ml dr^ rk Hr * + H dO^ dO^ + dz ^ dz 


(2.4a) 


where k is the coefficient of thermal diffusivity and has been non-dimensionalized with 
respect to c p . This is related to the Prandtl number a by 

(2.46) 

The dissipation function $ in (2.4a) is given by 
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(2.4c) 


2 K r dO ' dr r' J ' 3 dr 

The boundary conditions for the flow are that the no-slip conditions must be satisfied at 

the surface of the disc and there should be no motion far away from the disc. Hence, we 

have 

v r = v 9 = v z = 0 at z = 0, (2.5a) 


and 


v r — ► 0 and v$ -* — r as z — > oo. (2.56) 

The density, temperature, and viscosities must approach their free-stream values as z — *• oo, 
which, since they are non-dimensionalized with respect to these values, are equal to unity. 
We consider two possibilities for the boundary condition on the temperature at the surface 
of the disc. Either, the surface is adiabatic (i.e. thermally insulated) and then 


— — = 0 at z = 0, 
dz 


(2.6a) 
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or the surface is isothermal (i.e. the temperature of the surface is prescribed). In this case 
we suppose that the surface temperature is Tbw then the boundary condition is 

T = Tbw at z = 0. (2.66) 


3. THE BASIC FLOW 

We now consider the three-dimensional solution of equations (2.2) - (2.5) for the steady 
compressible boundary layer flow over a rotating disc, for which the Reynolds number will 
be large. We define the non-dimensional coordinates 

(r,0,z) = {r,0,R~$z B ), (3.1a) 


and velocities 

(v r ,v f ,v g ) = V B = (ub,v b ,R-*wb), (3.16) 

with the pressure, density, shear viscosity, bulk viscosity, conductivity, and temperature 
denoted by Pb>Pb> Mb>Mb>^b>T b , respectively. This flow will be defined by the governing 
equations (2.2) - (2.5) with d/dO = 0, since the flow will be axisymmetric. We substitute 
the basic flow variables (3.1) into the governing equations and neglect terms of 0(i2 -1 ) in 
the Navier-Stokes equations and 0(R“a) in the energy equation. Hence, for large Reynolds 
number the basic flow quantities are determined from the following equations: 


( du B , dt v Bs , dps , dp B Pb u b n 

+ ^ + UB -ai + + ~r °’ 


Pb(ub^- + wb 


oz B r 


+ JLfa**), 

dr dz B ' B dz B ' 


+ ^ + 2u fl ) = 

lr =0 ’ 

oz B 


Pb = 


(7 - mi 


d Tb , 
dz B 
1 d 


djpB 

dr 

dT B 


PbTb , 


" (..£ + + «k£>* + £?>*] 


dpB 

dz B 


dz B 
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(7 - 1 )M£ dz B dz B 
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(3.2a, 6, c,d,e, f) 


The boundary conditions are 


u B = v B = wb = 0 at z B = 0, 


(3.3a) 
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u B -*• 0, v B 


(3.36) 

(3.3c) 


— >• — r as z B — > oo, 

Pb,T b ,/x b -* 1 ,p B -^-pr as z B -* oo. 

Since we require some properties of the basic flow in the later stability analysis we proceed 
to find the solution of (3.2) and (3.3). For an incompressible flow the basic flow can be 
given by K&rmdn’s steady solution (see Kdrm&n’s (1921)). 

It is found that the present problem can be reduced to the corresponding incompressible 
one by use of the Dorodnitsyn - Howarth transformation which will be described below. 
(Riley (1964) obtains the same results using the Von Mises transformation.) Since the 
fluid is stationary everywhere outside the boundary layer the pressure gradient dp B /dr 
will be zero everywhere. Thus, from (3.2d,e) and the boundary conditions p B is constant 
and equal to ('•yM^,) -1 . We assume the Chapman viscosity law 

H = CT, (3.4a) 


where C is a constant, then the energy equation is decoupled from the momentum equa- 
tions. We take C = fx B w/T B w , where the subscripts w indicate values at the wall, then 
the law (3.4a) will be most realistic close to the wall. We use Sutherland’s law to evaluate 
fx B w > hence 

... ( i j .n. /r™ \ 

(3.46) 


_ T s / 2 I 1 + Cl / 1 " 00 \ 
“ BW \T B „ + Ci/Too j ’ 


where the constant C x is given by C\ = 198.6 °R. Since the pressure is constant from 
(3.2e) p B T B is constant and so the Chapman viscosity law (3.4a) gives p B fx B is constant. 
From the boundary conditions as z B — > oo we find that p B fx B = C. It is as a result of this 
expression that the energy equation is decoupled from the momentum equations. Thus, 
we have 

MB = {%=) T„ (3.5a) 


and 


Pb = 


l B 


(3.56) 


We make the Dorodnitsyn - Howarth transformation (see for example Stewartson 
(1964)) and introduce the coordinate 


1 r*B 

n =C&Jo 


Pb<Iz b , 


(3.6a) 


and the stream function ip, where 


1 dip 1 dip 

Ub = — x- and iv B = — , 

r op rp B or 


(3.66, c) 
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where d/ dr denotes a derivative with ip held constant. Then the equations (3.2) become 

ldip d , 13 ^. 1 dip d ld^\ v b _ & /1#V\ 

r dr} dr' r dr} r dr dr}'r drj r v & r dr} 2 'r dr} ’ 

1 dip dvB 1 dip dvs 1 dip 2 dip _ d 2 VB 

r dr} dr r dr dr} + r 2 drj VB + r dr} dr} 2 ’ 

1 _ 1 dip dT B s tfJLlHW )) 2 4 . (—B .\ 2 1 + 1 

(7 - 1 )Ml 0 r dr} dr r dr dr} "drj'rdr} 'dr] ' 0(7 — l JM 2 , cty 2 

(3.7a, 6,c) 


Now if we define ip by 

<P = -jH( 1 ,), 


then from (3.6b, c) 

= “#'(»?), (3.8a) 

C 1 / 2 

Wb = #(*?)» (3.86) 

Pb 

where a prime denotes differentiation with respect to rf. We also define the function G(rj) 
by 


v g = rG(r]). (3.8c) 

Hence, using (3.8) from (3.7a, b) we find that the functions H and G satisfy 

H"> + i^i! _ HH « _ 2 (G + l) 2 = 0, (3.9a) 

G" + H'[G + 1) - HG' = 0, (3.96) 

H'{ 0) = G(0) = 0, (3.9c) 

tf'(oo) = 0, G( 00 ) = -1. (3.9d) 


These are the equations satisfied by Kdrmdn’s steady solution which may be solved nu- 
merically. Two results from the solutions of (3.9) that we require in the later analysis 
are 

— H n (Q) 

u Q = = 0.51023 and v 0 = G'( 0) = -0.61592. (3.9c) 

z 

The energy equation (3.7c) becomes 


i- r 4-^r + ht 'b) = + 'W + 


n 


B 


(7-I )M 2 ) V 2 dr ' 4 v “ ; ' • > ' a (<y_i )M 2 3 ’ 

The solution to this equation can be found in terms of a heat conducting term and a 
viscous dissipation term. (See, for example, Riley (1964)). The solution can be written as 

Tb = 1 + ( T B w — l)^i (*?) + [l - 1) — 7 — <M*?)> (3.10) 
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where the functions fa and fa satisfy 


fa - oHfa = 0, (3.11) 

fa - oHfa + oH'fa = -o{{H") 2 + 4(G') 2 ). (3.12) 

For the special case of a model fluid with a — 1, exact solutions of fa can be obtained. 
Asymptotic solutions are given by Riley (1964) for large and small values of o. Note that 
if a is large the thermal boundary layer is much thinner than the viscous boundary layer 
and conversely, if a is small the thermal boundary layer is much thicker than the viscous 
boundary layer. In the analysis to follow we will consider a — 0.72, which is the value for 
air. For an adiabatic plate the boundary conditions are 

fa(0) = fa{0) = 0 , 

<M°o) = fa{oo) = o. 

Then, we find that the solution of (3.11) is fa = 0. For o — 1, the solution of (3.12) is 
given by 

h = 2 - ^ - 2(G) 3 , 

and if a — 0.72, numerically we find that fa(0) = 1.7899. While, if there is heat transfer 
at the wall then the boundary conditions are 


fa(0) = l,fa(0)=0, 

fa(oo) = fa(oo) = 0 . 


In this case, the solution of (3.11) is given by 


roo 

Mv) e~ cI dr), 

J " » 

( 3.13a) 

where 


IM = - f HWiri, 

Jo 

(3.136) 

A~ x = f°° e~ aI dr). 
Jo 

(3.13c) 


For o — 0.72, A is found numerically to be A = 0.3286. The solution of (3.12) has to be 
obtained numerically, except when o — 1 which gives 


* = -2(^ + G(G+ 1 )), 

and for a = 0.72 we find numerically that fa (0) = 0.9348. We will require the results for 
fa (0) , <t >2 (0) an d A in the later analysis. 
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4. THE WEAKLY NONLINEAR STABILITY ANALYSIS 


We consider the lower branch stationary viscous modes described by Hall (1986) and 
follow the analysis of MacKerrell (1987). We perturb the basic flow described in Section 3. 
The disturbance will have a triple-deck structure identical to that described by Hall (1986) 
for incompressible flow over a rotating disc. This will be based on the small parameter e, 
which is defined by 

£ = R-TS. 

The upper, main, and lower decks axe of thickness of 0(e 4 ),0(e 8 ) and 0(e 9 ), respectively. 
The structure of the disturbances in the main and upper decks is essentially the same as 
that found by Smith (1979), who investigated lower-branch disturbances to Blasius flow. 
It is found that the viscous effects are only important in the lower deck. The upper deck 
flow structure is essentially inviscid and irrotational and provides a pressure gradient which 
drives the lower deck. In the main deck the flow structure is inviscid with no variation in 
pressure across the layer. The lower deck is a thin viscous layer which is required to satisfy 
the no-slip boundary conditions at the disc surface. The disturbance is found to have 
wavenumbers in the radial and azimuthal directions of 0(e 4 ). Thus, we seek stationary 
! modes of instability with three-dimensional solutions proportional to 

E = exp(-^[f a(r)dr + 0O]). 

The wavenumbers a and (3 expand as 

0! = <*0 + G? OC\ + £ 3 Q!2 + • • * , 

! 

i /? = Po + + e s /?2 + • • • > 

i as in the corresponding incompressible problem. Hence, we add disturbances proportional 

, to E to the basic flow described in Section 3. The forms of the expansions of the disturbance 

will be different in each deck and are essentially fixed by balancing the convection and 
^ viscous terms in the radial momentum equation in the lower deck. 

We define the perturbed velocities, pressure, density and temperature by 

(v r ,v e ,v z )=V B + (U,V,W), (4.1a) 

{p,P>T) = ( Pb,Pb,T b ) + ( P,p p ,T p ), (4.1&) 

where the perturbation quantities are assumed to be small. The perturbation terms will 
| be in the form of fundamental terms and harmonic terms arising from the nonlineax nature 

| of the governing equations. Initially, we are seeking the solution at the position of neutral 

\ 9 

i 


t 


stability so a and /? will be real. Later we will consider the solution close to the position 
of neutral stability and then a will be complex. We begin by obtaining the solution in the 
upper deck. 


4.1. The Upper-Deck Solutions 

We introduce the coordinate Z where z = e 4 Z so that Z is 0(1) in the upper deck. The 
radial velocity perturbation expands as 

U = 6(e z U o + £ 4 U i + • • • )E + 6 2 {U 20 + cU 21 + • • ‘)E 2 + 6 3 ( — I- H — + • • •) J? 3 


+« ! (^r + + •••)« + o(« 4 ) + 


C.C. 


(4.2a) 

where c.c denotes complex conjugate of the preceding expression and 6 is a small amplitude 
which we shall later relate to e. The functions Uo,U 2 o etc. depend on r and Z. The 
azimuthal velocity, density, and temperature perturbations have similar forms to (4.2a). 
Note that there are no mean flow terms, i.e., terms independent of rand#. 

The vertical velocity and pressure perturbations have mean flow terms and are given 
by 

W = 6{e 3 Wo + e 4 W l + ''-)E + 6\W 2 o + eWn + ---)E 2 + 6 2 {e 5 W m o + e 6 W ml + -‘>) 
+<S (^ + ^ i + -) £S + « S (^ 2 + ^ i + -)® + 0(«*) + C .c. 

(4-26) 

P = 6{e 3 Po + s 4 P 1 + “‘)E + 6 2 {P 20 + eP n + --')E 2 + 6 2 {e«P m0 + e 7 7 ml + ••.) 

+ s * + ir + ' ' ') £S + (ir + it + -) B + °( {< ) + c - c - 

(4.2c) 

Using the Chapman viscosity law (3.4) the shear viscosity, bulk viscosity and thermal 
conductivity perturbations have the following form: 


Vbw rr 1 
Pv — rp T p> n p — 


p BWrp , _ 1 VBW rn 

ip, K p — i p . 

o -Ibw 


(4. 2d) 


From the solution for the basic flow described in Section 3, in the upper deck ub — 0,v B = 
—r,pB = 1,Tb = 1 and p B is constant. We substitute the perturbed flow (4.1) with the 
perturbed terms given by (4.2) into the governing equations (2.2)-(2.5). We first consider 
the linear solutions and equate terms proportional to 6E in the resulting equations. It is 
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found that the solution for the first order terms, which decays to zero as Z 

Fo = £/.-». r 0 = 4 - r 0 _ £/.-» 

7 0 = fe~ TZ , Po = Mlfe~ rz , T 0 = (7 - l)Mlfe~ TZ . 
Here f(r) is an amplitude function to be determined and 


00 , is given 


(4.3) 


r 2 



(4.4) 


where M r is a local Mach number defined by M r = rM^. Hence, from (4.4), for three- 
dimensional modes to exist for neutral wavenumbers we must have 


+ 4(1 - m?) > 0. 


(4.5) 


Thus, for \M r \ < 1 there will exist three-dimensional solutions for all real values of ao and 
/?o. However, for M r > 1 three-dimensional modes will only exist for T 2 > 0. We shall see 
later that from Hall (1986), in order for the three-dimensional stationary modes to exist 
the effective wail shear must vanish. Thus, we choose to look for stationary modes where 
the effective wall shear vanishes at zeroth order. The effective wall shear at zeroth order 
is 


du B , Podv B n A 

a 0 - 1 — = 0 at z B = 0. 

az B r az B 

Hence, from the basic solution described in Section 3 the condition (4.6a) is 


(4.6a) 


aor 

Ao 


« 1.207, 
«o 


(4.66) 


where Uo and Vq are defined by (3.9e) and the solution is neutral at the position r where 
(4.6) holds. This is the same result found by Hall (1986) for the incompressible problem. 
The conditions (4.5) and (4.6) imply that for the stationary three-dimensional modes of 
instability we are seeking to exist 

\M r \ < 1.567. (4.7) 


Thus, these modes will not be present in high Mach number flows. 

It is found that the next order terms in e have the same form as those given by (4.3) but 
with amplitude g(r). This unknown amplitude is determined from matching the expansions 
in the upper deck with the expansions in the main deck and will be a linear function of 
/(r). Analytic solutions can also be found for the higher order terms, but are not required 
in the present analysis. 
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If we now consider the terms proportional to 8 2 E 2 in the governing equations we find 
that the first harmonic terms are 

77 _ a ° f „-2TZ rr _h.-iYZ nr _ ** r .-2rz 

U20 — -7rj2 e j >'20— — e > W 20 — ~Z~h e > 

Po r Po 

P, o = he- 1TZ , = Mlhe-™, = (l - 1 )Mlf 2 e-* rz . 

Here / 2 (r) is an unknown amplitude function, which will later be found in terms of /(r). 
There is no contribution from the nonlinear terms until the solution for P 21 . It is found 
later that the harmonic terms proportional to 8 3 E 3 are not involved in the present analysis 
so we shall omit the solutions for these terms. Thus, we consider next the 0(6 3 E) terms 
which arise from the nonlinear interaction of the 0 (8E~ X ) terms with the 0(<5 2 £' 2 ) terms 
and the 0(<5£') terms with the 0(<5 2 £ , °) terms. In a similar way to the solutions given above 
we find that 

77 _ a ° f 0 -yz tt _f*.-TZ nr _ ^h-rz 

Uio — -p-j 4 e , V io — —e , w 10 — -p—e , 

Pio = he ~ rZ , p 10 = Mlf 4 e~ TZ , T 10 = (7 - l)M^/ 4 e- r ^, 

where f 4 is an unknown function of f(r). We note that the solutions for the velocities 
and pressure have the same form as those obtained for the incompressible problem by 
MacKerrell (1987), but where T now depends on M r as well as on a 0 and /? 0 . 

The solutions in the main deck must match with those obtained in the upper deck as 

Z^O. 

4.2. The Main-Deck Solutions 


From section 3 the boundary layer variable is defined by z B = e 8 z, where z B is 0(1) 
in this layer. The non-dimensional radial velocity perturbation has the form 




+ =£ + ...)*• + l*(=£ + + • ••) J + 0«‘) + c.c. 


(4.8a) 


with a similar expression for V . The vertical velocity perturbation has the following 
expansion: 

W = £(e 3 u> 0 + e 4 wi H )E + 8 2 (w 20 + H )E 2 + 8 2 (e 6 w m0 + e 7 w m i H ) 


+< 3 (^ + + ■ ■ -)E 3 + ^ ^ + ...)£ + 0(0 + c.c. 


,S, W J0 
-3 


wn 
2 


(4.86) 
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with a similar expression for P. The density and temperature perturbations are similar to 
each other. The expansion for the density perturbation is given by 


ft. = + + + + ^ + + + ^ + 


C2/P20 . />21 


2 ( PmQ Pml 


6 s 


+6*C^r + hjL + . . .)E 3 + s 3 ( ^ + hi + ...)J5 + 0(5 4 ) + c.c. (4.8c) 

We substitute the perturbed flow with the expansions (4.8) into the governing equations. 
Equating coefficients of the linear terms proportional to 6E gives the following solutions 
which match with the upper deck solutions as z B -* oo and Z — ► 0: 


Uq = 


r/dttB 

00 dz B ’ 


A 

Vo 


r f_dvB 

0o dz B ’ 


u>o 


*r / 

01 


(atQUg + 



Po = /> 


TfdRp 

0o dz B ’ 


r 

° 00 dz B ’ 


(4.9) 


where Ro denotes the density of the basic flow solution. The normal velocity wo satisfies 
the no-slip condition at z B = 0, but &o and Oo do not. As mentioned previously, from Hall 
(1986), in order to have stationary disturbances we must choose ao and 0o so that the wall 
shear is zero at zeroth order. Then aouo + 0oVo/r —*■ 0 as z B -*■ 0. At the next order of e 
the solutions which match with the upper deck solutions have the same form as (4.9) but 
with amplitude / replaced by g. Thus, we have pi = g. 

We find that the first order terms proportional to 6 2 E 2 have the same forms as (4.9) 
but with olq,0o and / replaced by 2olq,20q and / 2 , respectively. Similarly, the first order 
terms proportional to 6 Z E which match with the solutions in the upper deck have the same 
form as (4.9) but with the amplitude / replaced by / 4 . In both cases the nonlinear effects 
have no contribution until the third order solutions. 

We now require a lower deck to reduce the slip velocities to zero at the surface of the 
disc. The viscous terms in the governing equations will now play a part in the first order 
solutions at 0(5 E),0(6 2 E 2 ) and 0(5 S E I ) and the nonlinear effects will arise at the solution 
of the next order terms in e at 0 (6 2 E 2 ) and 0(5 3 2£). 
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4.3. The Lower-Deck Expansions 


We introduce the vertical coordinate £ where z = e g £ so that £ is 0(1) in the lower 
deck. For small values of zb we can expand the basic flow quantities u B ,v B ,Ro Then, in 
terms of £ 


u B = s£u B i + e 2 £ 2 «B 2 + 


v B = e£v B i + e 2 £ 2 t>B 2 + 


Ro = Row d* s £Roi d* £ 2 £ 2 -R<)2 d* ■ * ’ > T B = T B w + e£T B i + e 2 £ 2 T b 2 + • 
where Row = Rq( z B = 0). The other coeflicients are defined by 

1 d*u B I 


(4.10) 


U B j = 


i 1 




and similarly for v B j,Roj,T B j,j = 1,2, The coeflicients will depend on r and then the 

condition of zero wall shear can be written as a 0 u B i + PoV B i/r = 0. Hence, in order for 
the expansions of the disturbance in the lower deck to match with those in the main deck 
as £ — ► 00 , from (4.9) and (4.10) the radial velocity disturbance must have the form 

U = 6(— -^{u B \ + 2e£u B 2 + • • •) H — - + U\ + eU \ + • • ‘)E 

EPO £ 


+ 6 2 

+ 6 3 


((-Eh- 

\WVBl 


(«B1 + 2e£u B 2 H ) + ~^ ' Un 


H — jr + 
£■’ e 3 


■) 


£!+ CU + C^, + 


■) 




+0(6 4 ) *f- c.c. 

(4.11a) 

The azimuthal velocity, density and temperature perturbations have similar expressions to 
(4.11a) in the lower deck. The normal velocity perturbation has the form 

W = 6( — [{ a 0 U B2 + ~V B 2)£ 2 d" e(<XoU B s + — t>B3)£ 3 + ••*] + 6*Wq + E^Wi + •» -)E 

Po r r 


d*6 2 (( — '-Jjr — \{<XoU B 2 + —v B 2)i 2 H ] + e 3 W2 0 + e 4 W 2 i H )E 2 + s^WmQ + e 5 W ml H ) 

PO r 

+^ 3 (( — — + W 31 + '‘-)E 3 + ( ^J^[(o!ottB2 d" d" ■ ’ ’] 

E Epo r 


+W 10 + sW n + • • • )E ) + 0(6 4 ) + c.c. 


while the pressure perturbation is given by 


(4.116) 


I 
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P = 6{e s f + e 4 g + "')E + 6 2 ({f 2 + eP n + -")E 2 + e 10 P m o + e n P m i + ---) 

-|-6 3 ((_J£ + -j- H )E 3 + + — j- H )E) + 0(6 4 ) + c.c. (4.11c) 

c c c c 

Note that the mean flow terms in (4.11b) are larger than for the incompressible case. The 
perturbed flow given by (4.11) with the basic flow quantities given by (4.10) is substituted 
into the governing equations (2.2) - (2.5). 


4.4. The Linear Eigenrelation 


We equate coefficients of SE in the governing equations for the perturbed flow in the 
lower deck. It is found that 

-r/« fll U( 0 ,V 2 Al£) 

u ° - 5ST’ {412a) 

where U is a parabolic cylinder functions (see, for example, Abramowitz and Stegun (1964)) 
and 

« /?« r/V/t 

(4.126) 


a * , , 0 o ^ . ra 0 

A = 0 ! + ~ B! = 


from the r and 0 momentum equations of the basic flow solution and using the condition 
(4.6a) on the effective wall shear. Vq is determined from the following relation obtained 
from the continuity equation: 

oiqUo + —Vo = 0. 

r 

From the energy equation we find that 

-TfT Bl U{0,y/2eiAU) 

110 ~ ft U( 0,0) ’ 

and hence from the equation of state (2.3a) 

-r/iEoitt(0,-s/2a*Ai£) 


Pio 


01 


(4.12c) 


>o U( 0,0) 

We note that for an adiabatic wall Tio and p\o will be equal to zero. 

It is found that the functions Ui and V x cannot be determined until the solution for W 0 
is obtained. This is achieved by solving a third order ordinary differential equation similar 
to equation (4.16) in Hall (1986). Hence, Wo is given by 


Wo = — i[a.\UB\ + — vbi)~^t + A * 

r Po 


(rJl _|_ \ 

(«o + -j) 




Pbw 


+ 2*'« + ^) 


/?o\ tTSubM*) 


r 2 0 qTbwU(O,O)T B w 


ZiaorT f RqiF 3 {s) *(1 - a)T fRoirap F^s) 


ftpB\vU(0,0) 


20oPbwU(O,O) 


+ 


tT<jrra 0 ^ 2 
2^o TbwPbw ’ 


(4.12d) 
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where 5 = A* £. The function F\ satisfies the equation 

F"' - s 2 F[ + 2 sF 1 = 1, 

while F 2 ,F$ and F± satisfy the same equation as F 1 but with the right-hand-side replaced 
by U (0, V^2s) , £ [s^ (0, y/2<j "J s) ] and s 4 ^(0, \/2cf*s), respectively. The boundary conditions 
are 7^(0) = Fi(oo) = 0, i = 1 , 2 , 3 , 4 . The continuity equation relating U 2 and V 2 is 

./ tt . P o T/A . TT . Pi T/ \ , iTf , Pi . dW 0 ipioroco? . . 

t(cc 0 U 2 H V 2 ) + i(ctiU 0 -\ Vo) + H v.bi) + —77 = 0. (4.13) 

r r p 0 r d£ *Pbw 

Satisfying the boundary conditions at £ = 0 in (4.13) and using (4.12) gives the eigenrela- 
tion 

fn-2 I P°UF'(n\ 1 o-Yn-2 i rT f u Bi*l{0) Zia 0 rT f RoiF^O) 

(«o + )/*i(°) + 2t ( a o + -7) 7.3^ — 777—r 4,21,/n ^ 


r* J p*T BW U{ 0,0) 


PlU{ 0,0) 


t(l-a)r/i2oia 0 ri^(0) Tf{a 0 r)l 0 l 

— mam — = (1 + + t Bi *' 


(4.14) 


Since the right-hand-side of (4.14) is proportional to (l + t) and ao and Po are real we have 

ru B1 ri^(0) 


, P^F-tn) (1 ^r«„oorl!J(0) 

(“» + f) F - (0) - (1 - a) —, mrn~ [ ° + f) m 

ZoiorT RqiF$ ( 0 ) 


PgT BW U(0,0) 


P 2 oU(0,0) 

The values of F[( 0) and 273(0) /M (0,0) evaluated by Hall (1986) are given by 

7^(0) = 0.5991, 

22 ^( 0 ) 


(4.15) 


tt (0,0) 


= 0.4570. 


The values for 7^(0) and 7^(0) depend on a and can be obtained in the same way as Hall 
(1986). For a — 0.72 we obtain 

7^(0) = 0.0192, 

7^(0) = 1.6972. 

We shall consider solutions of the eigenrelation (4.15) for the cases of an adiabatic wall 
and an isothermal wall separately. 


4.4.1. An Adiabatic Wall 


In this case 72 qi = 0 so that the eigenrelation (4.15) becomes 




(4.16) 
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using (4.4). Note that if M r — 0 and Tbw = 1 the eigenrelation (4.16) reduces to the 
eigenrelation obtained by Hall (1986) for incompressible flow. Now let 


r 0 = 



(4.17a) 


From the solution of the basic flow described in Section 3, for an adiabatic wall the tem- 
perature at the wall is given by 


Tbw 


= 1 + (l ~ 1 ) 


MM 0) 


(4.176) 


where <fo(0) = 1.7899 for o = 0.72 and the value of 7 for air is 1.4. We also have from 
(3.6b) and (3.8) that 


ru 0 

UB1 ~ 7WCV2’ 


(4.17c) 


_ TV 0 

VB1 ~ T bw C V2’ 

where uq and Vo are given by (3.9e). The eigenrelation 
and (4.17) to give 


(4.17d) 

(4.16) may be rearranged using (4.6) 


r 0 r* = 


2u q F'(0) 

C'/>T* W FI(0)U(0,0) 




(4.18) 


Hence, from (4.18) Tors is a function of M r only. From the results obtained by Hall 
(1986) Tor^ = 1.224 for an incompressible flow where To can be regarded as the effective 
wavenumber of the perturbed flow. 

If we now equate real parts in the eigenrelation (4.14) we have, for real values of a 0 
and (3 0 , 

( a o + = 201tIw ( c * iUbi + ( 4 - 19 ) 

which can be rearranged using (4.17) to give 


&i G0P1 x * 

mr )r ‘ 


2(r 

|u 0 Uo|* 



(4.20) 


The value of (< x\/0o — otofii/ 0o)r* was found to be 2.312 by Hall (1986) for incompressible 
flow. The expression (4.20) gives a correction to the angle between the radius vector 
and the normal to the stationary vortices. The results given by (4.18) and (4.20) are 
the compressible generalization of the incompressible results obtained by Hall (1986). A 
discussion of these results is given in Section 5 for Too = 491.6°i2. 
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4.4.2. An Isothermal Wall 


If we manipulate (4.15) we obtain 


V 0 Hyf 2\ 4 f 


r 0 r* = 


RoiVoFi{0) , 2¥ 0 ^(0) 


(1 | p 0\ | SttpiZoiF^O) .!. , ■ 

2u 0 l/(0, 0) ' u§ ; Uolt(0,0) J 1 lUJ 


u, 


_1 

3. 


(4.21) 


From the solution for the basic flow 


M? 


2bi = [— (2bw ~ 1)-^ + (qf — 1)— p^2(0)]/(C 1 / 2 r jB jv), 

4 


where, for a = 0.72,44 = 0.3286 and <f> 2 (0) = 0.9348. From the equation of state Rqi = 
—TBi/Tg W . In order to evaluate (4.21) for values of M r the temperature of the wall must be 
specified. Tbw is the non-dimensional temperature at the wall, defined by Tbw = Tbw/T 
Hence, Tbw is the ratio of the dimensional temperature at the wall to the free-stream 
dimensional temperature. Calculations are obtained for Too = 491.6 °R and for a range of 
values of Tbw • Figure 1 shows Tor? as a function of M r for a — 0.72 for Tbw = 1.8, 1.0, 0.6 
and 0.2, and also for the case of an adiabatic wall. We see that the effective wavenumber 
of the disturbed flow decreases as the local Mach number increases. 

Equating real parts of (4.14) gives the same expression as for the adiabatic wall, i.e., 
(4.20), but of course the values of ap and /? 0 will be different, To being given by (4.21). 
Figure 2 shows (cti/Po — aoPi/Po)r* as a function of M r for a = 0.72 for Tbw = l-8> 1.0, 0.6, 
and 0.2, as well as for an adiabatic wall. We describe the orientation of the stationary 
vortices are given by the angle 0 O , where 0 Q is the angle between the tangent to the vortices 
and the normal to the radius vector, then 




tan(— - *o) 


ar 


apr 

Po 


+ e 2 r * 


A ( Ot\ &oPl\ & 

4 U pi ) r4 


+ 


Thus, from Figure 2 we see that d 0 decreases as the local Mach number increases for 
Tbw > 0-6 but for highly cooled walls with Tbw < 0.6 0 O increases as the local Mach 
number increases. The significance of these results is discussed in Section 5. 


4.5. The Nonlinear Solution 


We now consider the terms arising from the nonlinear interactions in the governing 
equations in the lower deck. We substitute the expansions (4.11) for the perturbed flow 
into the governing equations (2.2) - (2.5) and equate coefficients of 6 2 E 2 . We find the 


solutions 


r/ 2 umtt(0,v/2(2A)*Q 

2/?o U(0,0) 
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Ton — — 


r /tTbi U (0, y/2<r* (2 A)* £) 


z/(o,o) 


where V 2 o and p 20 are determined from 


aoU 20 + — V 20 = 0, 

r 

R01T20 — TbiP 20 = 0 > 

respectively. In the same way as for the linear terms, U 21 and V 21 cannot be determined 
until W 20 is obtained. The solution for W 20 is the same as that obtained for Wq, (4.12d), 
with a,-, f replaced by 2a,-, 2 / 2 , respectively, for j = 0, 1, 2. Satisfying the boundary 
conditions at the wall in the continuity equation gives the following expression relating / 2 
to /: 

, ,,,, , . . PL rru fll F;(0) _ S.aorrjio.F^O) 

W 4 (“0+ rI ) MB , + PIT BW U(0,0)V. BW «u sw U(0,0) 

1 r> P /3 . p2 /2 /3_ 


(4.22) 

Note that / 2 is proportional to f 2 . F$( 0) and F[( 0) may be eliminated from (4.22) by losing 
the linear eigenrelation (4.15). 

Now if the same procedure is followed for the solution of the terms proportional to 6 s E 
in the perturbation quantities solutions can be obtained for U 10 , Vio, W 10 , 2\o an d Pio- These 
are found to be identical to the linear solutions (4.12) with / replaced by the amplitude 
/t. However, the resulting eigenrelation here is 


k *r 7 2 


2ffi(M BW U{0,0) 


(l + .)r(q.r) 1 /» 

f^BW 1 20qTbw 


PL ru B1 rfj(0) 


(aiUB1 + S.„) + (.J + ^( 0 ) + 2 .(o| + 


3.TJi 01 a„r^(0) , i(l - «)r«,,a„rF|(0) ) .T »/*/,,. .. , A. , 

«o,o) + viuftof / = -^- (aoUB! + T” £!) ' 

(4.23) 

where * denotes complex conjugate. Since, from (4.22) / 2 is proportional to / 2 and from 
the linear eigenrelation (4.14) the left-hand-side of (4.23) is zero, this implies that / = 0. 
In order to determine whether the nonlinear effects are stabilizing or destabilizing we must 
move a sufficient distance from the position of neutral stability. 

We wish to determine an amplitude equation for the disturbance amplitude. This will 
involve derivatives with respect to r and nonlinear terms. In order to obtain the required 
equation, from the results of MacKerrell (1987) we find it is necessary to perturb the 
solution a distance of 0(e) from the neutral position. We must also choose 6 = e» so that 
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the appropriate nonlinear terms will arise in the eigenrelation. Thus, we define a new 
variable r\ by 

r = r + er u (4.24) 


where r is the position of neutral stability. Hence, we can write 

d _ ioco u*i 1 d 
dr e 4 e 2 + s dri ’ 

for the terms proportional to E and similarly for the 0(E 2 ) and 0(J? 3 ) terms. Thus, a 2 
may be effectively replaced by —id/dri in the equations previously obtained. The required 
eigenrelation will involve terms proportional to f,rif,df/dri and /|/| 2 . 

We replace r from (4.24) in the governing equations and rewrite the expansions of the 
perturbation quantities in the lower deck. For example, 

U = £*{^(ubi + 2e£v,B2 + • * •) + % + Ui + • • ’E 

+e3(J^(u.Bi + 2e£t*B2 4- • * *)*+ — ~ — H • • *)i? 2 +e2 (Umo + £ U m i + •••) 

+s(t^( u bx + 2e^«B2 H ) + H ) E 3 } + 0( ff4 ) + c * c * 

6Spo 

W = e»{( — — ^2 — [(ao u B 2 + —vb 2 )£ 2 + •••]+ e 6 Wo + s 7 W\ + •• -)E 
Po 5 r 

+ ^ ( _£M[ (aottfl2 + Po VB2 )f + ...] + e *w 2Q + . . . ) E 2 + e%{W m o + eW ml + • • •) 
Po^ r 

+e( — — — [(c*oUB2 + ~ u B2 )£ 2 + •••] + e 6 w s0 + • • ‘)E S } + 0( £l °) + c,c * 

Po r 

In order to obtain the required eigenrelation we need to determine W\. The linear terms 
may be obtained directly from the equation for W\. The terms occurring from expanding 
r are obtained from the 0(e) terms obtained by expanding r from (4.24) in the equation 
for W 0 . The nonlinear terms will simply be the nonlinear terms occurring in the equation 
for Wio. Then satisfying the boundary conditions at £ = 0 gives the solvability condition 


dW x 

di 


= 0 at £ = 0. 
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This condition leads to the following relation: 

■TV/ i q.rC‘/ ; f im . (ay + !3l)F[{0)^MlTl w C^ 
dr i 2fUoTfilnBW 1 \ r 2 u 0 (aor)2r 4 tZbr 3 

_ (1 , 3 . ^o 2 (agr 2 +^)^( 0 )Ti ty C 1 /2 (agr 2 + 

(aorj^uor (a 0 r)5F% 0 ri2oi 

(1 + 

(aor)Mo-RoiZ/(0,0) j 


, (1 + »KJ&rMBir0 1 ^ r~ 
(a 0 r)^r¥ 0 r 


s£/(0, 

^ (0, 0) 


££•£(£ = A-*s)ds. 


(4.25) 


Here 6 = (^ - l)M r 2 ^' 2 (0) / (4 T| h ,C 1 / 2 ) and is equal to zero for an adiabatic wall. The terms 
included in the expression RHS are linear terms depending on the functions Plo, Ali> Wo> U 0 , U u f 
and g and so can be written as 


i2HS(£)=c 1 (0/(r). 


We now follow the incompressible work of MacKerrell (1987) and write the eigenrelation 
(4.25) in the form 

= {C\r + iCu)f\f\* + ^l/(C2r + *@2%) + /(^Sr + *^3«)> (4-26) 

where C* r and C*, are constants for k — 1,2,3. Now multiply (4.26) by f* and add the 
complex conjugate to give 

^“(l/l 2 ) - 2C'ir|/| 4 + 2C 2r (ri + £~)l/| 2 * ( 4 - 27 ) 

In order to simplify (4.27) for C 2r ^ Owe introduce a new variable r[ = ri + C 3r / C 2r . Then 
(4.27) becomes 

Jr(l/r) = 2C 1 .|/| 4 + 2C 2r r' 1 |/| ! . (4.28) 

From (4.25) and the expression for / 2 , (4.22), after some manipulation we obtain 


where 

|A| 2 


C lr =f 
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3J^(o)(r*r 0 )i «5 _ m 2 ) s , v «| 4 (1 

2v / 2ui,r| l „|AP^ B , 1 r) wo 1 ’ 


I7(fir„)*(fi(0)) ! , |¥„tf„|fir„ 


st 2 


u 0 


«r2 


P’BW 


+ 


P'BW^'b 


BW 


(1 + + § - A4’)[(^i - ^)r-3] 2 

C v ug' u§ /?o /?0 
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lo^roj^fojluowol*/, i h , v l M 2\?( ai fta °w! 
V2T* bw S bw C'/> [1 + ul ) [1+ ul Mr) l /?0 * jr • 

Figure 3 shows 7 13 / 4 C lr as a function M r for an adiabatic wall and for an isothermal wall 
with T B w = 1.8, 1.0, 0.6, 0.2. We find that C\ r is always positive and thus the nonlinear 
effects are destabilizing. We find that 7 13 / 4 Ci r is less than the incompressible value of 
r 13 / 4 Ci r = 2.357 obtained by Mackerell (1987) for both cases of an adiabatic wall and an 
isothermal wall with Tbw > 1.0. However, for an isothermal wall with T B w < 1.0 r 13 / 4 C lr 
is greater than the incompressible value at M r = 0.05. Note that C lr decreases as the local 
Mach number increases for all values of Tbw • Thus, for a highly cooled wall the destabi- 
lizing effects of nonlinearity are greater for a compressible flow than for an incompressible 
flow. From (4.25) 


mn w (i + f >i(i + 1 - 

|u 0 u 0 |» «0 u o 


(r^r 0 )a(7- l)M r V a (0)ttg L * 

N*[(ZW ~ 1 )A - ( 7 - l)M r V' 2 (0)/4] V <) 



In Figure 4 r 7 / 4 C 2r is plotted as a function of M r for an adiabatic wall and for an 
isothermal wall with T B w = 1.8, 1.0, 0.6, 0.2, 0.15. We see that for the range of local Mach 
number where the solutions are valid and for T B w > 0.15 that C 2r is negative. Thus, in 
this case the amplitude of the disturbance increases or decreases depending on whether r 
is less than or greater than the neutral value. However, for T B w < 0.15 we find that C 2r 
becomes positive as M r increases. Thus, in this situation the amplitude of the disturbance 
increases as the distance from the neutral radius is increased. For an adiabatic wall and for 
an isothermal wall with Tbw > 0.6 we find that C 2r decreases as the local Mach number 
increases. However, for an isothermal wall with Tbw < 0.6 there will be a range of the 
local Mach number where C7 2r increases as the local Mach number increases. We note that 
the incompressible value of r 7 / 4 C 2r found by MacKerrell (1987) is equal to —0.903 and we 
see that for an adiabatic wall and an isothermal wall with Tbw > 1 that here r 7 / 4 C 2r is 
less than the incompressible value. 

With the substitutions 

|/| 2 = y|C 2r | 1/2 /(2C'i r )andri = x/\ 
forC 2r 7^ 0, (4.28) becomes 

^ - = ±2 xy + y 2 , (4.29) 
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with the positive sign taken if C 2r > 0 and the negative sign taken if C 2r < 0. If C 2r = 0 
use the substitutions 

|/| 2 = y (2C' 3r )/(2C7i r ) and = x/(2C Sr ), 


so that (4.27) becomes 



dy 2 

T = y + V- 
ax 

(4.30) 

The solutions of (4.29) and (4.30) are 



y{x) = e*7(l/|fo-/o* e^dt), C 2r > 0, 

(4.31a) 


y(z) = y 0 e x /(l + y 0 {l-e x )), C 2r = 0, 

(4.316) 


y(x) = e"*7(l/lfo “ (-A erfa: )/ 2 )> Cir <0, 

(4.31c) 

where 

yo = l/W = 0)l*15£ffa, C 2r ± 0, 

(4.32a) 

and 

So ■= \f(r[ = 0)P2C lr /2(CV), C Jr = 0. 

(4.326) 


The result (4.31c) is the generalization of the result obtained by MacKerrell (1987) where 
now C\ r and C 2r depend on Tbw and M r . Hence, for C 2r < 0 the same results hold for the 
compressible problem as for the incompressible problem. Thus, there exists a critical value 
of the initial amplitude of the disturbance. For values below this critical value (y 0 = 2 /a/tt) 
the disturbance will grow initially but will eventually decay to zero. Thus, we will have 
a stable solution. For initial amplitudes above this critical value the disturbance becomes 
exponentially large at a particular radius and does not decay. Thus, in this case the 
solution will be unstable. For C 2r > 0 we have growing solutions terminating at finite 
values of the radius. Thus there are no stable solutions for C 2r > 0. Figure 5 shows y as a 
function of x for different values of t/o from (4.31c) for C 2r < 0. 

From the results of MacKerrell (1987) for an incompressible flow r _19 / 8 (— C' 2r ) 1 / 2 /(2C'i r ) = 
0.202. Figure 6 shows r -19 / 8 |/(0)| 2 /y 0 = r" 19 / 8 (— Cj,.) 1 / 2 /^*?!,.), for C 2r < 0, as function 
of M r for an adiabatic wall and for an isothermal wall with Tbw = 1.8, 1.0, 0.6, 0.2. Thus, 
we see that r -19 / 8 (— C 2r ) 1 / 2 /(2C lr ) is greater than the incompressible value for an adiabatic 
wall and for an isothermal wall with Tbw >1.0 but less than this value for a range of the 
local Mach number when Tbw < 1.0. 

Thus, for an adiabatic wall and an isothermal wall with Tbw > 1.0 we require a 
larger value of the initial amplitude of the disturbance in the compressible case than in 
the incompressible case for the solution to become unstable. However, for some values 
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of the local Mach number for an isothermal wall with Tbw < 1.0 we require a smaller 
value of the initial amplitude of the disturbance in the compressible case than in the 
incompressible case for the solution to become unstable. Clearly, from Figure 6, as the 
local Mach number increases the threshold amplitude is increased and thus we require a 
greater initial amplitude for the disturbed flow to become unstable for larger values of the 
local Mach number. 

5. CONCLUSION 

The main conclusion to be drawn is that the effects of compressibility greatly alter 
these stationary viscous modes of instability since they cannot exist for large enough Mach 
number. We find that the results for a compressible flow, when these stationary modes 
exist, are similar to those obtained by MacKerrell (1987) for an incompressible flow. The 
main difference occurs for highly cooled walls. For instance, the finite amplitude of the 
disturbance solution increases or decreases depending on whether r is less than or greater 
than the neutral position. In the same way as for the incompressible problem for particular 
values of M r and T BW we have the idea of a threshold solution as described by Stuart (1971) 
for the finite disturbance amplitude. Thus, there exists a critical value of the amplitude 
of the initial disturbance below which the disturbance eventually decays to zero as the 
distance from the position of neutral stability is increased and the solution will be stable. 
However, if the initial amplitude of the disturbance is above the critical value the solution 
will end in a singularity at a finite radius and the solution will be unstable. For highly 
cooled walls ( Tbw < 0.15) this threshold solution does not exist thus we have only unstable 
modes of instability. The nonlinear effects are found to be destabilizing for a compressible 
flow, which is also the case for an incompressible flow (see MacKerrell (1987)). It is found 
that for an adiabatic wall and an isothermal wall with T B w >1.0 the effects of nonlinearity 
are not as important here as for an incompressible flow and decrease as the local Mach 
number increases. But for an isothermal wall with Tbw <1.0 the effects of nonlinearity 
can be more important for a compressible flow than for an incompressible flow. 

The weakly nonlinear analysis has enabled the effect of Mach number on the finite 
amplitude growth of a disturbance close to the neutral location to be described. Two cases 
of an adiabatic wall and an isothermal wall were investigated. It was found that the nature 
of the results was the same for an adiabatic wall and an isothermal wall with Tbw > 1.0. 
For these cases it was determined that the wavelength of the disturbance increases as the 
local Mach number increases. In addition, as the temperature of the wall is raised for an 
isothermal wall the disturbance wavelength increases. The orientation of the stationary 
vortices is described by the angle 0 O where 0 O is the angle between the tangent to the 
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vortices and the normal to the radius vector. As the local Mach number is increased we 
find that 0 O decreases. For an isothermal wall as the temperature of the wall is increased 0 o 
decreases. Thus, we observe that for an isothermal wall with Tbw >1.0 increasing either 
the local Mach number or the temperature of the wall produces the same effect. For an 
isothermal wall with Tbw <1.0 the same results apply with the exception that 0 Q may 
increase as the local Mach number increases. 

We now describe the main differences between our results for a compressible flow and 
those obtained by Hall (1986) and MacKerrell (1987) for an incompressible flow. We find 
that there is a larger mean flow correction in the normal direction in the lower deck for 
a compressible flow than for an incompressible flow, but otherwise the disturbance flow 
structures are the same. In Figure 1 we see that IV 1 / 2 is less than the value obtained by 
Hall (1986) for an adiabatic wall and also for an isothermal wall with Tbw > 1.0. Thus, the 
wavelength of the disturbance is longer in these situations than for an incompressible flow. 
However, for an isothermal wall with Tbw <1.0 the wavelength of the disturbance may be 
shorter than that for an incompressible flow. Figure 2 shows that (ai//? 0 — ao^i/^ 2 )r 5 / 4 
is greater than the incompressible value obtained by Hall (1986) for an adiabatic wall and 
an isothermal wall with Tbw > 1.0. Hence, the angle of orientation of the stationary 
vortices 0 O is smaller for these cases than for an incompressible flow. For an isothermal 
with Tbw < 1.0 0 O may be larger than for an incompressible flow. 

From the nonlinear analysis we found that r ls / 4 C lr is less than the value obtained by 
MacKerrell (1987) for an incompressible flow for an adiabatic wall and for an isother- 
mal wall with Tbw > 1.0. Thus, in these cases the nonlinear effects are less important 
for a compressible flow than for an incompressible one although they are still destablish- 
ing. However, for an isothermal wall with Tbw < 1.0 r 13 / 4 Ci r may be greater than the 
incompressible value in which case the nonlinear effects will be stronger than for an in- 
compressible flow. From Figure 4 we find that for highly cooled flows C 2r may be zero or 
positive. Thus, in this case there will be no threshold solution with only unstable modes of 
instability present. From Figure 6 we see that for C 2r < 0 r~ 19 / 8 \f(r[ = 0)| 2 /y 0 is greater 

than the value obtained by MacKerrell (1987) for an adiabatic wall and for an isothermal 
wall with Tbw < 1.0. The solution (4.31c) breaks down when y 0 = 2/y/n which is the 
same value as for an incompressible flow. Hence, for an adiabatic wall and an isothermal 
wall with Tbw > 1.0 for an unstable solution the initial amplitude of the disturbance for a 
compressible flow must be larger than that required for an incompressible flow. But for an 
isothermal wall with Tbw < 1.0 we require a smaller initial amplitude for the solution to 
become unstable. In addition, for a highly cooled wall there is the possibility of no thresh- 
old solution, the solution being unstable in this case. Thus, an unstable mode of instability 
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is more likely to occur for a compressible flow with a highly cooled wall than for an in- 
compressible flow. For an adiabatic wall or an isothermal wall with Tbw > 1.0 the system 
must be forced even more strongly for a compressible flow than for an incompressible flow 
for these stationary modes of instability to be evident. 

We found that these stationary viscous modes of instability do not occur above a critical 
supersonic local Mach number of about 1.57. This suggests that for small amplitude 
disturbances for values of the local Mach number greater than 1.57 the inviscid stationary 
modes described by Gregory, Stuart and Walker (1955) will be more dominant. Hence, 
for small amplitude disturbances it is unlikely that this mode would be observed in flows 
over bodies moving at very high Mach numbers. On the other hand, for values of the local 
Mach number less than 1.57, for highly cooled walls there is a strong possibility that the 
nonlinear effects may cause the early breakdown of the laminar flow. 
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Figure 1 . The effective wavenumber for an adiabatic wall and 
for an isothermal wall with Tbw = 1 .8, 1 .0, 0.6, 0.2 and 
also showing the corresponding incompressible value. 



t 


I 

I 

I 

I 



5 

in 


cQlIcm o 

a"|ca? 


I 


Incompressible 


C\J 

o 




CQ 



CO 


CO 





o 


o 

"O C/3 

®T3 

> § 
S'n 

~ O 
5 CD 

.2 o 
*o - 
cc q 

c i - 1 

ro ,rC 

i_ co 
o • 


© II 
03 £ 

I? 

II 

■2 03 
c E 
.5 o 
o £ 
0) o 

C/) 


o 

o 


c 

03 


CD 

iSo 


cvi 

CD 

V- 

3 

g> 

Ll 


showing the corresponding incompressible result. 



Figure 3. r 13/4 C 1r for an adiabatic wall and for an isothermal wall with 
Tbw = 1 .8, 1 .0, 0.6, 0.4 and also showing the corresponding 
incompressible result. 


Tbw = 0-15 



Figure 4. r 7/4 C 2 r as a function of Mr for an adiabatic wall and for an 
isothermal wall with T B w = 1 8, 1 .0, 0.6, 0.2, 0.15 and also 
showing the corresponding incompressible result. 



Figure 5. The behavior of y for a range of values of y 0 . 
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Figure 6. The initial disturbance amplitude for C2r < 0 for an adiabatic wall 
and for an isothermal wall with Tew = 1 .8, 1 .0, 0.6, 0.2 and also 
showing the corresponding incompressible value. 
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